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Ever since Karl Schwarzschild’s 1916 discovery of the spacetime geometry describing the interior 
of a particular idealized general relativistic star — a static spherically symmetric blob of fluid 
with position-independent density — the general relativity community has continued to devote 
considerable time and energy to understanding the general-relativistic static perfect fluid sphere. 
Over the last 90 years a tangle of specific perfect fluid spheres has been discovered, with most of these 
specific examples seemingly independent from each other. To bring some order to this collection, 
in this article we develop several new transformation theorems that map perfect fluid spheres into 
perfect fluid spheres. These transformation theorems sometimes lead to unexpected connections 
between previously known perfect fluid spheres, sometimes lead to new previously unknown perfect 
fluid spheres, and in general can be used to develop a systematic way of classifying the set of all 
perfect fluid spheres. 
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I. INTRODUCTION 


General relativistic perfect fluid spheres, or more precisely general relativistic static perfect fluid spheres, are 
interesting because they are first approximations to any attempt at building a realistic model for a general relativistic 
star U, 2j, fe m- The central idea is to start solely with spherical symmetry, which implies that in orthonormal 
components the stress energy tensor takes the form 

- p o o o - 

rp _ 0 Pr 0 0 

0 0 p t 0 1 J 

.0 0 0 p t _ 


and then use the perfect fluid constraint p r = pt, making the radial pressure equal to the transverse pressure. By 
using the Einstein equations, plus spherical symmetry, the equality p r = pt for the pressures becomes the statement 

G§§ = Gr? = ^ 00 - ( 2 ) 

In terms of the metric components, this leads to an ordinary differential equation [ODE], which then constrains the 
spacetime geometry, for any perfect fluid sphere. 

Over the last 90 years, many “ad hoc” approaches to solving this differential equation have been explored, often by 
picking special coordinate systems, or making simple ansatze for one or other of the metric components |5i,jf|, 7, : 8j. (For 
recent overviews see ms .) The big change over the last several years has been the introduction of “algorithmic” 
techniques that permit one to generate large classes perfect fluid spheres in a purely mechanical way fgj, [jfj [ll| ■ In 
this article we will extend these algorithmic ideas, by proving several solution-generating theorems of varying levels 
of complexity. We shall then explore the formal properties of these solution-generating theorems (which generalize 
the notion of the Buchdahl transformation) and then will use these theorems to classify some of the previously known 
exact solutions, and additionally will generate several new previously unknown perfect fluid solutions. 
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II. SOLUTION GENERATING THEOREMS 

Start with some static spherically symmetric geometry in Schwarzschild (curvature) coordinates 

ds 2 = -C(r) 2 di 2 + + r 2 dO 2 (3) 

B(r) 

and assume it represents a perfect fluid sphere. That is, Ggg = G?? = G 11 . While Ggg = G^ is always fulfilled due 
to spherical symmetry, setting G?? = Ggg supplies us with an ODE 

[r(rQ']B' + [2r 2 C" - 2{rQ'}B + 2C = 0 , (4) 

which reduces the freedom to choose the two functions in equation © to one. This equation 0 is a first-order linear 
non-homogeneous equation in B(r). Thus — once you have chosen a ((r) — this equation can always be solved for 
B(r). Solving for B(r) in terms of £(r) is the basis of ESE1, (and is the basis for Theorem 1 below). On the other 
hand, we can also re-group this same equation as 

2r 2 BQ" + ( r 2 B' - 2rB)Cf + (rB' -2B + 2)( = 0, (5) 

which is a linear homogeneous second-order ODE for £(r), which will become the basis for Theorem 2 below. The 
question we are going to answer in this section is, how to systematically “deform” the geometry 0 while still retaining 
the perfect fluid property. That is, suppose we start with the specific geometry defined by 

ds 2 = -Co (r) 2 d t 2 + + r 2 dfl 2 (6) 

Bo(r) 

and assume it represents a perfect fluid sphere. We will show how to “deform” this solution by applying five different 
transformation theorems on {Co, .Bo}, such that the outcome still presents a perfect fluid sphere. The outcome of this 
process will depend on one or more free parameters, and so automatically generates an entire family of perfect fluid 
spheres of which the original starting point is only one member. In addition, we analyze what happens if you apply 
these theorems more than once, iterating them in various ways. 


A. Four theorems 


The first theorem we present is a variant of a result first explicitly published in El , though another limited variant 
of this result can also be found in Q. We first re-phrase the theorem of El in slightly different formalism, and 
demonstrate an independent way of proving it. Using our proof it is easy to show that by applying theorem 1 more 
than once no further solutions will be obtained, therefore the transformation in theorem 1 is, (in a certain sense to 
be made precise below), “idempotent”. 

Theorem 1. Suppose (Co( r ), Bp(r)} represents a perfect fluid sphere. Define 


A oM 


CoW 


Co(r) +r CoM. 



Co(r) Co(r) ~ r Co (r) , \ 
Co( r ) Co( r ) + r Co( r ) J 


Then for all X, the geometry defined by holding CoM fixed and setting 


( 7 ) 


ds 2 = -Co(r ) 2 d t z + 


dr 2 


B 0 (r) + A A 0 (r) 

is also a perfect fluid sphere. That is, the mapping 

7~i(A) : (Co, Bq} i— ► {Co, B 0 + AA 0 (Co)} 


dfr 


( 8 ) 

( 9 ) 


takes perfect fluid spheres into perfect fluid spheres. Furthermore a second application of the transformation does not 
yield new information, T 1 is “idempotent”, in the sense that 

Ti(X n ) o • • • o Ti(\f) o Tj(Ai) : {Co,Bp} >—» jCo,Bp + 1 Xj] A 0 (Co)| (10) 

We also note that Ti always has an inverse 


[T.(A )]- 1 = U- A). 


( 11 ) 
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Proof for Theorem 1. Assume that {Co(t)i -Bo(r)} is a solution for equation 0. Under what conditions does 
{(o(r), Bi(r)} also satisfy equation 0 ? Without loss of generality, we write 

B\{r) = B 0 (r) + A A 0 (r). (12) 

Equation 0 can now be used to determine Ao (r). That ordinary inhomogeneous first-order differential equation in 
B now simplifies to 


[r-K 0 )']A' + [2r 2 ^-2(rCo)'] A 0 = 0, 


(13) 


which is an ordinary homogeneous first-order differential equation in A 0 . A straightforward calculation, including an 
integration by parts, leads to 


A„(r) = 


[Ko)'] 


exp 


K 

Ko)' 


dr 


(14) 


Adding and subtracting ±2(rCo 2 )/(CoKo)') t° the argument under the integral leads to 


An = 


Co (r) 


Co(r) +r Co(r) 


2 

r exp 


l 1 Co(r) 


oW Co{r)~rCo(r) 


Co (r) Co(r)+rCoW 


dr 


(15) 


as advertised. If we apply this transformation a second time we obtain no additional information. To see this, consider 
the sequence 


{Co,B 0 }^{Co,B 1 }^{Co,B 2 }... (16) 

But at the second step, since Co has not changed, we have Ai(r) = Ao(r). More generally, at all subsequent steps, 
Aj(r) = Ao(r). We can write this as 

. (17) 

or in the more suggestive form 

n 

n T i=Ti (is) 

i=1 

where the symbol = indicates 11 equality up to relabelling of the parameters ”. That is, transformation Ti is “idempotent” 

up to relabelling of the parameters. (See figure 0) □ 

A version of Theorem 1 can also be found in Jjjj]. Specifically, after several manipulations, changes of notation, 
and a change of coordinate system, the transformation exhibited in equation (16.11) of [§) can be converted into the 
transformation Ti of Theorem 1 above. 

Applying theorem 1 to a fixed {CoKo} generates a one-dimensional space of perfect fluid spheres, which leads to 
the corollary below: 

Corollary 1. Let {Co,-Ba} and {Co,-B&} both represent perfect fluid spheres, then for all p 

{Co,pB a + (l-p)B b } (19) 

is also a perfect fluid sphere, furthermore all perfect fluid spheres for a fixed Co can be written in this form. 

Proof. The result is automatic once you note that for fixed Co the ODE for B is linear, (though not homogeneous, 
which is why the two coefficients p and 1 — p are chosen to add up to 1) . □ 

We defer extensive discussion of the application of this theorem and its corollary until section EH where we use 
this and our other generating theorems as a basis for classifying perfect fluid spheres. At this stage we mention, only 
as a very simple example, that Ti applied to Minkowski space results in the geometry of the Einstein static universe. 


n T i(A,) = T, hr 


i= 1 


K=1 
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FIG. 1: The solid lines show B(r ) for 10 reapplications of Theorem 1 onto the Minkowski metric. The dashed line corresponds 
to a single application with the specific choice A onC e = A;. It can be seen that 10 applications of Theorem 1 can be expressed 

by one application. 


Theorem 2. Let {Co,So} describe a perfect fluid sphere. Define 


Zo(r) = a + i 


r dr 


CoM 2 \ABoO) 

Then for all a and e, the geometry defined by holding Bo(r) fixed and setting 

ds 2 = -Co(r ) 2 Z 0 (r ) 2 df 2 + + rW 

B 0 {r) 

is also a perfect fluid sphere. That is, the mapping 

T 2 {&,c) ■ (Co, So} | —*■ {Co i?o(Co,So),5o} 


( 20 ) 


( 21 ) 


( 22 ) 


takes perfect fluid spheres into perfect fluid spheres. Furthermore a second application of the transformation does not 
yield new information, T 2 is “idempotent” in the sense that 


T 2 {a n ,e n ) o ■■■ o 72 ( 0 - 3 , 63 ) o 72 (cr 2 ,e 2 ) o 7 . 2 ( 0 - 1 , 61 ) = T 2 (a n . .. a 3 a 2 cr 1 , e n ... 32 i), 


(23) 


where 


6n...321 = (eio-20-3 ■ • ■ cr n ) + (0-1 1 e 2 a 3 ■ ■ ■ a n ) + {cr 1 1 a 2 1 e 3 • • ■ a n ) H-b (a ± 1 a 2 V 3 1 • • • e n ). (24) 

Furthermore, theorem 2 is invertible (as long as a ^ 0): 

[72(o-, e)]^ 1 = 7 2 ( 1 /ct, -e). 


(25) 
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Proof for Theorem 2. The proof of theorem 2 is based on the technique of “reduction in order”. Assuming that 
{( 0 (r), B 0 (r)} solves equation 0, write 


CiO) = CoM Z 0 (r). (26) 

and demand that (Ci(r), B 0 (r)} also solves equation 0. We find 

(r 2 toB' 0 + 4r 2 C 0 B 0 - 2r^B 0 )Z' 0 + (2r 2 Co B 0 )Zg = 0 , (27) 

which is an ordinary homogeneous second-order differential equation, depending only on Zq and Z'f. (So it can be 
viewed as a first-order homogeneous differential equation in Z ', which is solvable.) Separating the unknown variable 
to one side, 


M = 2 <o 1 

Z' 0 2 B 0 Co r‘ 

Re-write Z'f /Z'q = dln(Zb)/dt, and integrate twice over both sides of equation (12811 . to obtain 

_ rdr 

zj o — (7 6 - . 

J Co (r) 2 y/B 0 (r) 

depending on the old solution {Co(r), -Bo( r )}; and two arbitrary integration constants a and e. 
To see that the transformation T 2 defined in Theorem 2 is “idempotent” we first show 


T 2 o T 2 = T 2 , 


(28) 


(29) 


(30) 


and then iterate. The precise two-step composition rule is 


T 2 (ct 2 , e 2 ) 0 T 2 (cti, ei) =T 2 


(^ 201 , 


£102 


£2 

cri 


(31) 


To see “idempotence”, note that for fixed Bq{t) equation 0 has a solution space that is exactly two dimensional. 
Since the first application of T 2 takes any specific solution and maps it into the full two-dimensional solution space, 
any subsequent application of T 2 can do no more than move one around inside this two dimensional solution space 
— physically this corresponds to a relabelling of parameters describing the perfect fluid metric you are dealing with, 
not the generation of new solutions. To be more explicit about this note that at step one 


Co —Ci — Co { a i + £1 


dr 


while at the second step 


Cl —C 2 = Cl \ + £2 


C 0 W 2 \jBo(r) f 


r dr 

C 1 M 2 V B o (r) 


That is: 

C 2 = Co { a i + £1 


r dr 


dr 


Co (r) 2 y/B 0 (r) 
But this can be rewritten as 


< J 2 + £2 


Co(r) 2 y/B 0 (r) [01 + £i Jr dr/(Co(r) 2 yjB 0 {r)f 


£2 / dZo 


C 2 = Co ^0 { <72 + - I } = Co Z 0 <1 / d ( — ) } = Co Zo <j 0 2 - - 


T, Hi 


1 1 
Z 0 &1 


Therefore 


C 2 = Co <{ -- + 

£l 


, £2 1 

o 2 -\ - 

£1 <?\ 


Zo>- 


(32) 


(33) 


(34) 


(35) 


(36) 

















from which the composition law 


T 2 (a 2 ,e 2 ) ° T 2 (cti, ei) = T 2 


(38) 


^cr 2 cri, 


. e 2 

eicr 2 H- 

o-i 


follows immediately. (Note that the composition law for T 2 is actually a statement about applying reduction of order 
to second-order ODEs, it is not specifically a statement about perfect fluid spheres, though that is how we will use it 
in the current article.) The general composition law then follows by induction. To find the inverse transformation we 
choose ct 2 = 1 /ct i and ei = —e 2 , for which 


T 2 (1/<7 i, — ei) o T 2 (oi, ei) =T 2 ( 1 , 0 ) = I. 


Comment: As other special cases of the composition law we also mention the results that 

n / n \ 

ni 2 (Mo=T 2 i,5> , 


i— 1 


i— 1 


and 


\2{&,€) = I 2 l <7 5 £ <7 ~r CT • • • + (T v ’ + G K ’ I. 

Now as long as er > 1 then for sufficiently large n we see 

T 2 (ct, e) n » T 2 (<t", ^e) = a”" 1 T 2 (a, e) 4 T 2 (a, e), 


(39) 

□ 


(40) 


(41) 


(42) 


where at the last step we have used the fact that the overall multiplicative factor a n ~ 1 can simply be reabsorbed into 
a redefinition of the time coordinate. Because of this result, we see that for fixed a > 1 (and fixed but arbitrary e) 
repeated numerical applications of T 2 (cr, e) will have a well-defined limit. In figure[5]we have tested the composition 
law numerically. 

A strictly limited version of theorem 2, with little comment on its importance, can be found in [f|. Specifically, 
after several manipulations, changes of notation, and a change of coordinate system, the transformation exhibited in 
equation (16.12) of |S| can be seen to be equivalent to the sub-case a = 0, e = 1 of transformation T 2 above. 

For some purposes it is more useful to rephrase theorem 2 as below: 

Corollary 2. Let {£ a , Bq} and {0>, So} both represent perfect fluid spheres, then for all p and q 


{pCa + q Cb, B 0 } 


(43) 


is also a perfect fluid sphere. Furthermore, for fixed Bq all perfect fluid spheres can be written in this form. 

Proof. The result is automatic once you note that for fixed Bq the ODE for ( is second-order linear and homogeneous. 

□ 

We again defer extensive discussion of the application of this theorem and its corollary until section mu at this 
stage mentioning only as a very simple examples that T 2 applied to either the Einstein static universe or the anti-de 
Sitter universe results in Schwarzschild’s stellar solution (position-independent density). Similarly, corollary 2 applied 
to a combination of the Einstein static universe and anti-de Sitter space is another way of obtaining Schwarzschild’s 
stellar solution. 

Having now found the first and second generating theorems makes it useful to define two new theorems by 
composing them. Take a perfect fluid sphere solution {<(o,-Bo}. Applying theorem 1 onto it gives us a new perfect 
fluid sphere {^o, B i}. The new B\ is given in equation 0. If we now continue by applying theorem 2, again we get 
a new solution {<(,1?!}, where C now depends on the new B\. All together we can consider this as a single process, 
by introducing the following theorem: 






7 



FIG. 2: The solid lines show £(r) for 10 reapplications of Theorem 2 onto the Minkowski metric. The dashed line corresponds 
to a single application with the specific choice for a once and e onC e determined by the composition law of Theorem 2. It can be 
seen that 10 applications of Theorem 2 can be re-expressed by a single application. 


Theorem 3. //{Co, Bo} denotes a perfect fluid sphere, then for all cr, e, and X, the three-parameter geometry defined 
by 


ds 2 = -Co(r) 2 {cr + e 


r dr 


Co (r ) 2 \/B 0 {r) + A A 0 (r) 


dr + 


dr 2 


B 0 (r) + A A 0 (r) 


r 2 dfr 


(44) 


is also a perfect fluid sphere, where A 0 is 


A o (r) 


CoM 

Co(r) +r Co(r) 



Co<» Co{r)-rCo(r) 1 

Co(r) Co(r)+rCoW J ' 


(45) 


That is 


T 3 —T 20 T 1 : {Co, -B 0 } > {Co, B 0 + A Aq(Co)} >—> {Co ^o(Co,B 0 + AA 0 (Co)), B 0 + AA 0 (Co)}- (46) 


Now, instead of starting with theorem 1 we could first apply theorem 2 on {Co,-Bo}. This gives us a new perfect 
fluid sphere {Ci, B 0 }, where Ci = Co ^o(Co, B 0 ) is given by equation J2J). Continuing with theorem 1 leads to {Ci, B}, 
where B depends on the new Ci- Again, we can consider this as a single process, by introducing the following theorem: 

Theorem 4. //{Co, Bo} denotes a perfect fluid sphere, then for all a, e, and A, the three-parameter geometry defined 
by 


ds 2 


-Co(r) 2 



r dr 

CoM 2 \/Bo(r) 


dr -I_ 

™ + B 0 (r)+ A A 0 (Cr,r) 


r 2 dfl 2 


(47) 
























(48) 


is also a perfect fluid, sphere, where Ao(Ci,r) is defined as 

Ci(r) ~ r f'flr) 1 

Ci(r)+r C((r) J 

depending on = Co Zq, and where as before 

Zo(r) = a + e f - r -±= . 

J Co(r -) 2 V B o(r) 

That is 


Ao«i ,r) = 


Ci (r) 


Ci(r) + rC((r) 


r 2 exp 


CiM 


(49) 


T 4 — Ti o T 2 : {Co, B o} 1-1 • {Co ^o(Co, -So), -So} > {Co ^o(Co, -So), -So + AA 0 (Co ^o(Co, -So))}- (50) 

Theorem 3 and theorem 4 are in general distinct, which can be traced back to the fact that theorem 1 and theorem 
2 do not in general commute. Furthermore theorem 3 and theorem 4 are in general not idempotent 

T 3 #T 4 ; T 3 oT3#T 3 ; T 4 oT 4 #T 4 . (51) 

The best way to verify this is to try a few specific examples. There may be some specific and isolated special metrics 
for which theorem 3 and theorem 4 happen to be degenerate, or idempotent, and finding such metrics is important 
for our classification programme. 


B. Formal properties of the generating theorems 


The solution generating theorems we have developed interact in a number of interesting ways and exhibit numerous 
formal properties that will be useful when classifying generic perfect fluid spheres. To start with, theorem 3 (T 3 ) and 
theorem 4 (T 4 ) can be expressed in terms of theorem 1 (T 3 ) and theorem 2 (T2): 

T 3 ( 5 ):=(T 2 oT 1 ) 5 = T 2 (T 1 ( 5 )); 

T 4 (s):=(T 1 oT 2 ) 3 = T 1 (T 2 ( 5 )). 1 j 

where g is a metric representing a perfect fluid sphere. Having this in mind, some of the new solutions generated by 
starting from some specific solution can be identified. For example: 

T40T1 = TioT 2 oTi = Ti o T 3 , (53) 


or 


T30T3 = T 2 0T1 oT 2 0T1 = T 2 0T4 0T1. 


(54) 


The idempotence of T 3 and T 2 in this formalism is: 

(Ti°Ti)g=Ti(Ti(g)) = Tl ( 3 ); 
(T 2 oT 2 ) g =T 2 (T 2 (g)) = T 2 (g). 


(55) 


Taken together, it is possible to simplify all formulae wherever Ti and T 2 appear next to each other more than once. 
The following examples should demonstrate how this works: 

(T 2 O T 3 ) 3 = (T 2 o T 2 o Ti) «7 a (T 2 o Ti) 5 = T 3 ( 5 ); 

a (56) 

(Ti oT 4 )g= (Ti O Tx O T 2 ) g A ( Tl O T 2 ) g = T 4 ( 5 ), 


and in the same way 


T40T3 = Ti oT 2 oT 2 0T1 = Ti oT 2 0T1 = TioT 3 = T4oTi5 
T30T4 = T 2 o Ti o Ti o T 2 — T 2 oTi oT 2 = T 2 oT 4 = T 3 oT 2 . 


(57) 


These relationships can be used to structure the solution set generated starting from any particular perfect fluid sphere, 
and moreover to classify which metrics can be produced by our theorems, and ones which cannot. For example, the 
idempotence property of theorem 1 and theorem 2 enables us to divide the class of perfect fluid spheres into seed 
metrics and non-seed metrics. Seed metrics can never be generated by using one of the two theorems Ti or T 2 , while 
non-seed metrics are connected to other simpler metrics via one of these theorems. We formalize this in the following 
subsection. 
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C. Seed and non-seed metrics 

Definition (Seed metric): Take a metric g (or a parameterized class of metrics) and apply theorem 1 or theorem 
2 on it. Three different cases are possible: 

• Each of the applications supplies us with a new solution. [Ti (g) ^ g ^ T 2 (< 7 ).] We define a metric with this 
behaviour as a seed metric. (We shall soon see several examples of this behaviour.) 

• Only one of the applications supplies us with a new solution, while the other one gives us the same metric we 
started with. [Ti(c/) = g or T 2 (g) = g.) These metrics are non-seed metrics. (We shall soon see several examples 
of this behaviour.) 

• Both applications give us the same metric we started with. [T 1 (g) = g = T 2 ( g ).] Metrics of this type are fixed 
points of the transformation process and we then also have T 3 (g) = g = T^g). 

While we have encountered numerical examples that seem to exhibit this fixed point behaviour, we have no analytic 
proof for the existence of non-obvious fixed-point metrics. There is one obvious but not particularly useful example 
of a fixed point class of metrics. If we take the ODE in equation 0 , and write down its most general solution as 
a functional of the arbitrary parameters C( r )i then any of our solution generating theorems applied to this most 
general solution will at most move us around in the parameter space characterizing the most general solution - — the 
most general solution of equation 0 , or equivalently equation ©, is thus an infinite-parameter fixed point of the 
generating theorems. But apart from this obvious example, it is unclear whether other fixed point classes of metric 
exist. 

We can nevertheless develop several formal lemmata regarding fixed-point metrics. For instance 
Lemma 1 . Suppose we have a metric such that V cr, e, A 

Ts(o', e,X)g = g, (58) 

where we recall that = denotes equality up to redefinition of parameters. Then in particular 

T iW g - g - T 2 (a,e) g, (59) 

and so 

T 4 (a, e, X) g = g. (60) 

Conversely, suppose we have a metric such that 

T 4 {<j, e, X) g = g, (61) 

then 

Ti{\) g = g = T 2 (a,e) g, (62) 

and so 


Tjfoe.A) g = g- (63) 


Proof. Trivial, note that T 3 (cr, e, A = 0) = T 2 (ct, e) and T 3 (ct = 0, e = 0, A) = Ti(A). 

Then recall T 4 = Ti 0 T 2 . 

Similarly for the converse. 

□ 

Lemma 2. Suppose we have a metric g such that 


T 4 g — T 3 g 

(64) 

and then define g' by 


T 4 g — g = T 3 g. 

(65) 

Then we have 


14 g = g = T 3 g 

( 66 ) 

so that g' is a “fixed point” of both T 3 and T 4 . 



10 


Proof. Note 


T 4 g' 4 T 4 o T 3 g 4 T x o T 2 o T, g A Tl o T 3 g A T, o T 4 g A T 4 3 = </ 


(67) 


and similarly for T 3 : 


T 


3 3 


= T 3 o T 4 5 4 T 2 o T, o T 2 a A T 2 o T 4 g A T 2 o T 3 g A T 3 g = g' 


( 68 ) 

□ 


Several other formal theorems along these lines can be constructed, but these seem the most important results. 
Finally, among the formal properties enjoyed by the generating theorems, we mention the fact that theorems 3 and 
4 are “conjugate” to each other in the following sense 



—1 

III 

H 

0 

—1 

to 

— T 1 0T2 0 Ti 

°[Ti] 

7 

t 

0 

CO 

O 

1 — 

II 

1 

(69) 

and similarly (when the 

appropriate inverse 

[T2] — 1 exists) 





T 3 = T2 0 Ti 

= T2 0 Ti 0T2 

°[t 2 ]- 

1 

Osl 

0 

h- 

0 

CNI 

1 - 

II 

1 

(70) 

We can write this as 








T 3 (ct, e, A) ~ 

T 4 (ct, 

e, A), 

(71) 


where ~ denotes the concept of “similarity” under conjugation by invertible generating theorems. We wish to empha¬ 
sise that similarity ~ is a statement that holds for particular and fixed values of the parameters (er, e, A), as opposed 
to = which denotes equivalence under redefinition of parameters. If one is working numerically, it is much easier to 
ask questions involving similarity For analytic work, it is typically easier to ask questions involving equivalence =. 


D. Two linking theorems 


The last two solution generating theorems we shall present are slightly different from those developed so far: Consider 
a perfect fluid sphere solution {Co,-Bo} and extend it to a new perfect fluid sphere {(qZ 0 ,B 0 + A}. Previously, we 
had either set Z 0 = 1 and obtained theorem 1, or we had set A = 0 and obtained theorem 2. In other words, we only 
changed one metric component at a time. (From this point of view theorem 3 and theorem 4 are, strictly speaking, 
not new theorems, in that they are replaceable by iterations of theorem 1 and theorem 2 and vice versa.) We now 
investigate what happens if we place no a priori restrictions on Z and A, and allow both metric components to vary 
simultaneously. The differential equation 0 for this problem now becomes 

[r(rCo^o)'] A' + [2r 2 (Co^o)" - 2(rCo^o)'] A + [r 2 Co^ + dr 2 C 0 B 0 - 2r( 0 B Q ] Z' 0 + 2r 2 C Q B Q Z'f = 0. (72) 

Note that if A = 0 this becomes equation Q, while if Zq = 1 this becomes equation (|27|). In general, this ODE 
of first-order in A, and — as long Z is not a constant — inhomogeneous. In terms of A this ODE can be solved 
explicitly and the result stated as a new theorem: 

Theorem 5. Suppose {(jo,-Bo} describes a perfect fluid sphere, and let Zo(r) be arbitrary. Define 


A(A, r) 


AoM 


A- 


[(4r 2 qB 0 + r 2 Co B' 0 - 2rf 0 B 0 ) Z' Q + 2 r 2 C 0 B 0 Z"] (Co + r^} 2 

r3 (r( 0 Z 0 y Co 2 Z 2 


exp 



(Co^o)' Co^o - r ((qZq)' ,1,1 
Co^o Co^o + r (Co^o)' d J d J’ 


(73) 


where 


Ao (r) 


[ C 0 Z 0 V o f 0 f (W CqZq - r (CoZo)' , ) 
VCo^o + r iCoZo)') r XP \ 2 J CaZ 0 ( 0 Z 0 + r (Co^o)' ' 


Then for all A, the geometry defined by an arbitrary chosen Zo(r) and setting 


dr 2 


(74) 


ds 2 = -Co (r) 2 Z 0 (r) 2 dt 2 + 


Bo{r) + A(A, r) 


+ r 2 clD 2 


(75) 
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corresponds to a perfect fluid sphere. That is, the mapping 

75 (A) : {Co, -Bo} | — 1 ’ (Co Z Q , B 0 + A(A, Co)} (76) 

takes perfect fluid spheres into perfect fluid spheres. 

Note that if Zo(r) = 1 this simply reduces to theorem 1. Re-arranging equation Cl in terms of Zq leads to a 
second-order inhomogenous differential equation, which cannot in general be solved for a prescribed A, unless one 
imposes further constraints. So further exploration in that direction is moot. There is however a related theorem 
(which may be easier to understand) in terms of parametric derivatives: 

Theorem 6 . Let {C (p), B(p)} denote a one-parameter class of perfect fluid spheres , so that the differential equation 
0 is satisfied for all p. Then 

MO'! {%) + - 2 «» (f) + (|)" + ^ ^ (|)’ + W - 2B + 2 > (|) = °- < 77 > 

In particular if dt/dp = 0 this reduces to the ODE II.Hi . while if dB/dp = 0 this reduces to the ODE 0). 

This is simply an alternative viewpoint on the previous theorem, emphasising the differential equation to be solved. 
We again defer extensive discussion to the next section, but that this stage point out that if we invoke theorem 5 
and apply it to Minkowski space, then making the choice Zq = 1 + r 2 /a 2 leads to the general Tolman IV metric — 
that is: 

T 5 (Minkowski; Z 0 = 1 + r 2 /a 2 ) = (Tolman IV). (78) 

Even before we systematically start our classification efforts, it is clear that the solution generating theorems we have 
established will inter-relate many of the standard perfect fluid spheres. 


E. Formal properties of the linking theorems 


Before turning to issues of systematic classification of perfect fluid metrics, we wish to establish a few formal 
properties of the linking theorems. To simplify the notation, let us define the differential expression 

£>(C, B) = [r(rC)']B' + [2r 2 C" - 2{rC,)']B + 2( = 2r 2 Bf" + ( r 2 B' - 2 rB)f + (; rB ' -2 B + 2)C- (79) 

Then the condition for a perfect fluid sphere is simply 

D(C,B)=0. (80) 

Now define 

D((o,Bo + Ao) = D(Co, Bo) + -Di((o; Ao). (81) 

Then it is easy to check that 

Di((o-, Ao) = [r(rCo)'] A' + [2r 2 ^ - 2(rCo)'] A 0 . (82) 

The ODE for theorem 1 , where we assume {Co,-Bo} is a perfect fluid sphere, is then 

-D 1 (Co; A 0 ) = 0. (83) 

Now define 

D(toZo , -Bq) = Z 0 D(to , -Bo) + D 2 {f 0 , B 0 ; Z Q ), (84) 


then it is easy to check that 

D 2 (Co, -Bq; Z 0 ) = (r 2 t 0 B' 0 + 4r 2 (’ 0 B 0 - 2r( 0 B 0 )Z' 0 + (2r 2 Co-B 0 )^o ■ (85) 

The ODE for theorem 2, where we assume {Co,-Bo} is a perfect fluid sphere, is then 


D 2 (Co, -Bo; Zq) — 0. 


(86) 
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Now let us consider any simultaneous shift in ( and B 1 as considered in theorem 5. We have 

-D(Coi?o, B 0 + A 0 ) = D((oZo, B 0 ) + Di((oZo; A 0 ) = Z 0 D((o, B 0 ) + D 2 (Co, Bo', Z 0 ) + Di(CoZ 0 ; A 0 ). (87) 

But now let us write 

-Di(Co-£o; A 0 ) = ZqD^q- A 0 ) + D 12 (Co; Z 0 , Ao), (88) 

where a brief computation yields 

^ 2 (Co; Z 0 , A 0 ) = r 2 Co^A' + [2r 2 Co^' + 4r 2 C^ - 2rCo^]A 0 . (89) 

Then all in all 

■D(Co^Oi B 0 + Ao) = Z 0 D(( 0 , B 0 ) + Z 0 Di((q; A 0 ) + -Di 2 (Co; Zo, A 0 ) + -D 2 (Co, B 0 ; Z 0 ). (90) 

So if Co^ Bq and (Co-Zo, -Bo + Ao} are both perfect fluid spheres we must have 

^o£>i(Co; Ao) + -Di 2 (Co;^o, A 0 ) + D 2 (( 0 , Bq; Z 0 ) = 0. (91) 

This is 


Z 0 x (ODE for theoreml) + (cross term) + (ODE for theorem2) = 0. (92) 

The cross term vanishes if either Zq = constant or Aq = 0 in which case we recover the usual theorem 1 and theorem 
2. If we do things in the opposite order, then 

D((oZo, Bo + Ao) = ZqD((o,Bo + Ao) + D 2 { Co, Bq + Ao; Zq); (93) 


— Z 0 D(( 0 , Bo) + ^oDi(Co, A 0 ) + -D 2 (Co, B 0 + A 0 ; Z 0 ). (94) 

We now have to compute 

B 2 ((o,B 0 + A 0 ; Z 0 ) = D 2 (Co: B 0 ; Z 0 ) + D 2 i(Co, B 0 ; Zq, A 0 ), (95) 

and after a brief calculation 

D 2 i(Co, -So; Z 0 , Ao) = 2r 2 CoA 0 i?o + [^CoAq + 4r 2 CoA 0 — 2rCoA 0 ]^o = Di 2 (Co, B 0 \ Z 0 , Ao). (96) 

Thus the cross term is the same, no matter how you calculate it, and we still have the identity 

D(Co^o, B 0 + A 0 ) = Z 0 D(( 0 , Bo) + Z 0 Di((o', A 0 ) + Di 2 (Co; Zo, A 0 ) + D 2 ( Co, B 0 ;Z 0 ). (97) 

Provided Co, B 0 and {Co-Zo, B 0 + A 0 } are both perfect fluid spheres we again deduce 

^oDi(Co; A 0 ) + Di 2 (Co; Zq, A 0 ) + D 2 (Co, -B 0 ; Zo) = 0. (98) 

Now this gives us another way of looking at theorem 3 and theorem 4. For theorem 3 we first apply theorem 1 so we 
have the two equations 


Di(Co; A 0 ) — 0; 


and 


D 2 (Co, B 0 + A 0 ; Zo) = -Di 2 (Co; Z 0 , A 0 ) + D 2 (( 0 ,Bq; Z 0 ) — 0. 

Conversely, for theorem 4 where we first apply theorem 2 we see that we need to solve 

D 2 (Co,B 0 ; Z 0 ) = 0; 

and 


(99) 

( 100 ) 

( 101 ) 


Di(Co^o; Ao) = Z 0 D i(Co; A 0 ) + Di 2 (Co; Z 0l Ao) — 0. (102) 

For theorem 5 we pick Zq arbitrarily, and solve the single ODE 

iMCo^o; Ao) + D 2 (Co, Bo] Z 0 ) = 0. (103) 

This is a single first-order linear inhomogeneous ODE for Ao, and hence solvable. (In particular this makes it clear 
that theorem 5 is an inhomogeneous version of theorem 1 with a carefully arranged “source term” Z) 2 (Co, Bq', Zo). 
While theorem 5 is not “idempotent” it does satisfy the important formal property that: 
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Lemma 3. 


T: o T 5 A t 5 , (104) 

which in particular tells us that the output from theorem 5 is never a seed metric. 

Proof. Applying theorem 5 we need to solve 


Di{CqZ 0 \ A 0 ) + .£> 2 (Co, B 0 ; Z 0 ) — 0, (105) 

in order to map 

{Co, -B 0 } —> {Co- 2 ' 0 , B 0 + A 0 }. (106) 

Now apply theorem 1 to (Co-^o, -Bo + Ao}, this means we have to solve the homogeneous ODE 

D 1 (Co^o;A 1 ) = 0. (107) 

But then, using properties of first-order ODEs 

UMCo^o; Ao) + D 2 (Cq, Bq; Z 0 ) = 0} © {Di(Co-Zo; Ai) = 0} (108) 

=> {P>i((oZ 0 ; Ao + Ai) + D 2 ((o,Bo; Zq) = 0 } (109) 

which is the ODE from theorem 5 back again. □ 


(The net result of this observation, as we shall see in the next section, is that theorem 5 can be used to connect 
one seed metric with the “descendants” of another seed metric.) 

III. CLASSIFYING PERFECT FLUID SPHERES 

We will now see the power of these transformation theorems (solution generating theorems) by using them in a 
number of different ways: to generate several new perfect fluid spheres, to relate various perfect previously known fluid 
spheres to each other, and to classify the geometries we encounter. First some minor comments and warnings: Despite 
comments made in [lj, Kuch 2 1= Tolman V; and Kuchb I b is a perfect fluid for general values of its parameters. 
Furthermore RR-I = Einstein static; RR-V = Tolman V (n = —5/4). If we had not noted these degeneracies then 
our tables below would have been more extensive, but would have conveyed no extra information. 

Starting with the metric for any known perfect fluid sphere and successively applying theorem 1 and theorem 2 
numerous times, will supply us with endless “new” perfect fluid sphere solutions. Some of these “new” solutions 
might already be findable in the literature, some of them might be truly novel. Some of these solutions can be 
written down in fully explicit form. Some solutions are explicit but not elementary, in the sense that while the metric 
components can be exhibited as specific and explicit integrals, these integrals cannot be done in elementary form. 
Some solutions are so complex that present day symbolic manipulation programs quickly bog down. (For specific 
symbolic computations we have used a vanilla installation of Maple.) 

To summarize the situation we present several tables and diagrams. Two tables are used to provide the names and 
explicit metrics for the perfect fluid spheres we consider. Two other tables are used to describe the inter-relationships 
of these perfect fluid spheres under Ti, T 2 , T3, and T4. In these tables the notation “[integral]” means that it is 
definitely a novel perfect fluid solution but the metric components involve an explicit integral that does not appear 
to be do-able by elementary methods, and is so complicated that it does not seem worthwhile to even write it down. 
Recall that a seed metric is one for which theorem 1 and theorem 2 both yield metrics distinct from the seed: 
Ti(ff) ^ g $= T 2 (<7). In contrast for non-seed metrics one or the other of these theorems is trivial, either Ti(g) = g or 
T 2(g) = 9 . 

By considering the idempotence of theorem 1 and theorem 2, and the fact that theorem 3 and theorem 4 can be 
expressed in terms of the first two theorems, it is possible to structure, and therefore to graphically visualize the 
relationship between all metrics generated from a given seed metric. We demonstrate this behaviour starting with 
Minkowski spacetime as seed metric, checking if it is possible to create endless new perfect fluid spheres starting with 
this trivial seed. The first few steps can be carried out explicitly, and show that Minkowski space generates several 
well known interesting perfect fluid models. 
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Some selected perfect fluid solutions 


Name 

Metric 

Minkowski 

Einstein static 

de Sitter 

Schwarzchild Interior 

Schwarzchild Exterior 

—dt 2 + dr 2 + r 2 AQ 2 

-At 2 + (1 - r 2 /I? 2 ) -1 dr 2 + r 2 d0 2 

— (1 — r 2 /R 2 ) dt 2 + (1 — r 2 /R 2 )- 1 dr 2 + r 2 dfl 2 

-{A- B^l-r 2 /R 2 ^) dt 2 + (1 - r 2 /7? 2 ) -1 dr 2 + r 2 dfi 2 

— (1 — 2m/r) 2 dt 2 + (1 — 2m/r) 1 dr 2 + r 2 dQ 2 

SI 

K-O III 

Kuchl lb 

B-L 

—r 4 dt 2 + dr 2 + r 2 dD 2 

— {A + Br 2 ) 2 dt 2 + dr 2 + r 2 dD 2 

— (Ar + Br In r) 2 dt 2 + 2d r 2 + r 2 dfl 2 
-A(r 2 /a 2 )dt 2 + 2(1 + r 2 /a 2 )~ 1 dr 2 + r 2 dfi 2 

Tolman IV 

Tolman IV (b —> oo) 

~B 2 ( 1 + r 2 /a 2 ) dt 2 + dr 2 

B 2 (1 + r 2 /a 2 ) dt 2 + ^^dr 2 + r 2 dfl 2 

Tolman V 

Tolman V (A -> 0) 
Tolman VI 

Tolman VII 

Tolman VIII 

B 2 r 2 (l+n ) dt 2 + (2 n 2 ) [1 j4r 2(2-nO/(2+n)]-l dr 2 + ^^2 

_ B 2 r 2(i+n) dt 2 + (2 _ n 2 )dr 2 + r 2 dD 2 

-( Ar 1-71 + Br 1+n ) 2 dt 2 + (2 - n 2 )dr 2 + r 2 dfl 2 

— B 2 cos |ln ^1 — 2r 2 /a 2 + r 4 /b 4 + r 2 /a 2 — b 2 /a 2 J + #| 

+(1 — 2 r 2 /a 2 + r 4 /6 4 ) -1 dr 2 + r 2 dfl 2 

/[ 2 r 2 (n-l)(n- 4 )/n ( -n 2 , Br ~(n 2 -8n+8)/n . Q r ~ 2 (n 2 - 4 n+ 2 )/n\ 

l (n 2 — 4 n+ 2 )(n^ — 8 n+8) ' ' J 

+ ((„*-4„ + 2 ) (^- 8 „ +8) + Br-^-M/r, + Cr —2(„ 2 —4n+2)/n^ - 1 ^2 + ^2 

Kuch 68 II 

Kuch 68 I 

M-W III 

Kuch I b 

Heint Ila (C=0) 

Heint Ila 

— (1 — 2 m/r) dt 2 - 

— ^Ay/l — 2 m/r - 

+(1 — 2 m/r)~ 
Ar(r a) + 

— r 2 [A + B arctan] 

— (1 + ar 2 ) 3 dt 2 + 

— (1 + ar 2 ) 3 dt 2 + 

- [(1 - 2 m/r) (1 + C(2r - 2m) 2 )] 1 dr 2 + r 2 dQ 2 

- B |^r 2 /m 2 + 5 r/m — 30 + 15i/l — 2m It ln{l — r/m — i Jr(r — 2m)/m}j ^ dt 2 
" 4 dr 2 + r 2 dfl 2 

sy^ydr 2 + rW 

l (a/y/a 2 + r 2 )] 2 dt 2 + 2(1 + r 2 /a 2 )^ 1 dr 2 + r 2 dfl 2 
(1 + ar 2 ) /(I - ar 2 /2)]dr 2 + r 2 dQ. 2 

1 WTT—Z) + Cr /-1 dr 2 + r 2 dfl 2 

^ i +“ r J (l+a7-2)Vl+4 ar 2 


TABLE I: Some well-known perfect fluid spheres and their corresponding metrics. Note that we have often reparameterized 
these metrics to make them easier to deal with, and so their appearance (but not the substance) may differ from other 
sources Qii Metric SI is a special case of K-O III, Tolman V, and Tolman VI notable for its extreme simplicity. 


In figure 0 all the ideas from the previous sections are used. Each box represents a specific metric (perfect fluid 
sphere) while the arrows correspond to the application of the different theorems. The horizontal arrows correspond to 
an application of theorem 1 and the vertical arrows to an application of theorem 2. In addition, the vectors pointing 
along the diagonals can either indicate an application of theorem 3 or theorem 4. A dotted arrow corresponds to the 
application of theorem 3 while the dashed arrow represents an application of theorem 4. 

Figure [3] shows that after applying theorem 1 to the Minkowski seed metric , we get the Einstein static universe. By 
the idempotence of theorem 1, n applications of Ti to the seed metric still results in the Einstein static. Similary, any 
number of applications of theorem 2 after the first (which leads to the K-0 III solution) does not give us any further 
new solutions (see figure 0. We can also see that the first and second generating theorems are not commutative. 
Application of theorem 1 and theorem 2 in that order to the Minkowski seed metric results in the Schwarzchild Interior 
geometry, whereas application of theorem 1 after theorem 2 gives us the Martin 3 solution, and the three-parameter 
Schwarzschild interior solution is quite distinct from the three-parameter Martin 3 solution. 

Indeed, it seems as if it is possible to create endless new solutions for a perfect fluid sphere out of the Minkwoski 
metric (or any other of the seed metrics). After several iterations the calculations become more complex, and can 
no longer be carried out analytically. We then resort to numeric computation to find out whether theorem 3 and 
theorem 4 have some sort of numerical limit, a numerical fixed point, or not. Depending on the choices made for 
the parameters (A, cr, e) both theorems converge very quickly. Specifically, we used Matlab to numerically analyze 
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Some apparently new perfect fluid solutions 



TABLE II: Some apparently new perfect fluid spheres and their corresponding metrics. Sometimes the relevant integrals cannot 
be done in elementary form. When they can be done they are explicitly shown. 

Some seed geometries and their descendants 


Seed 

Theorem 1 

Theorem2 

Theorem 3 

Theorem 4 

Minkowski 

exterior Schwarzscliild 

de Sitter 

Tolman V (A = 0) 

SI 

M-W III 

Heint Ila (C=0) 

Einstein static 

Kuch68 II 

Tolman IV 

Tolman V 

Tolman V (n = +1) 
Martin 2 

Heint Ila 

K-0 III 

Kuch 86 I 

interior Schwarzscliild 

Tolman VI 

K-O III 

P5 

P3 

interior Schwarzschild 
[integral] 

P7 

Wyman III 

P4 

P6 

P8 

Martin 3 
[integral] 
[integral] 
Wyman Ila 

Martin 3 
[integral] 
[integral] 


TABLE III: Seed solutions and their generalizations derived via theorems 1-4. The notation “[integral]” denotes a metric so 
complicated that explicitly writing out the relevant integral is so tedious that it does not seem worthwhile. 


the evolution of £(r) and B(r) after applying theorem 3 several times to the Minkowski metric. Numerical results 
are summarized in figure 01 This figure indicates that theorem 3 and theorem 4 both appear to have well defined 
numerical limits, though we have no direct analytic solution for the fixed point metric. 
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Some non-seed perfect fluid geometries and their descendants 


Base 

Theorem 1 

Theorem2 

Theorem 3 

Theorem 4 

Tolman IV 

B-L 

Heint Ha 

Tolman IV 

B-L 

Heint Ila 

P7 

Kuchb I b 

P8 

P7 

Kuchb I b 

P8 

[integral] 

[integral] 

[integral] 

Tolman VI 

Kuchl lb 

K-O III 

Wyman Ila 

Martin 1 

Martin 3 

Tolman VI 

Kuchl lb 

K-O III 

[integral] 

P2 

PI 

[integral] 
Martin 1 

Martin 3 


TABLE IV: Non-seed solutions and their generalizations. 



FIG. 3: Structure graph for Minkowski space as seed metric. 

IV. DISCUSSION 

Using Schwarzschild coordinates we have developed several transformation theorems that map perfect fluid spheres 
into perfect fluid spheres, and have used these transformations as a basis for classifying different types of perfect fluid 
sphere solutions. If we apply these theorems on a known perfect fluid sphere, different solutions are often obtained. 
While some of the solutions are already known in the literature, most of them are novel. Therefore, we have developed 
a tool to generate new solutions for a perfect fluid sphere, which does not require us to directly solve the Einstein 
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FIG. 4: Numerical implementations of Theorem 3 and Theorem 4 can be used to investigate their behaviour for n applications, 
as the number of applications becomes large. For specific choices of parameters A;, at and ti successive applications of Theorem 
3 and Theorem 4 appear to be approaching a limit, in the sense that both metric components seem to converge toward fixed 
points. 


equations. (Of course the whole procedure was set up in such a way that we are implicitly and indirectly solving 
the Einstein equations, but the utility of the transformation theorems is that one does not have to go back to first 
principles for each new calculation.) 

In addition, we have also established several relationships among the generating theorems. Previously, all metrics 
seemed to have nothing more in common than representing a perfect fluid sphere. We have developed the concept of 
a seed metric, which is one that cannot be generated by our theorems. Starting with a seed metric and applying our 
theorems, it is possible to structure, therefore to visualize, the relationship between all metrics generated from a given 
seed metric in a graph. Based on this example it is possible to create endless new solutions out of the Minkowski 
metric. We also used a numeric program to investigate whether theorem 3 and theorem 4 have fixed point limits. 
Both seem to converge very quickly. 

In summary, the class of static perfect fluid sphere solutions to the Einstein equations exhibits a deep and somewhat 
unexpected structure — the transformation theorems of this article provide a new and distinct way of looking at the 
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problem, and yield a new way of viewing the inter-relationships between different static fluid spheres. 


[1] M. S. R. Delgaty and K. Lake, “Physical acceptability of isolated, static, spherically symmetric, perfect fluid solutions of 
Einstein’s equations,” Comput. Phys. Commun. 115 (1998) 395 arXiv:gr-qc/9809013 . 

[2] M. R. Finch and J. E. F. Skea, “A review of the relativistic static fluid sphere”, 1998, unpublished. 
http://www.dft.if.uerj.br/usuarios/JimSkea/papers/pfrev.ps 

[3] H. Stephani, D. Kramer, M. MacCallum, C. Hoenselaers, and E. Herlt, Exact solutions of Einstein’s field equations , 
(Cambridge University Press, 2003). 

[4] D. Martin and M. Visser, “Bounds on the interior geometry and pressure profile of static fluid spheres,” Class. Quant. 
Grav. 20 (2003) 3699 arXiv:gr-qc/0306038 . 

[5] H. A. Buchdahl, “General Relativistic fluid spheres”, Phys. Rev. 116 (1959) 1027-1034; 

“General Relativistic fluid spheres II: general inequalities for regular spheres”, Ap. J. 146 (1966) 275-281. 

[6] H. Bondi, “Spherically Symmetrical Models In General Relativity”, Mon. Not. Roy. Astron. Soc. 107 (1947) 410. 
“Massive spheres in general relativity”, Mon. Not. Roy. Astron. Soc. 282 (1964) 303-317. 

[7] M. Wyman, “Radially symmetric distributions of matter”, Phys. Rev. 75 (1949) 1930. 

[8] S. Berger, R. Hojman and J. Santamarina, “General Exact Solutions Of Einstein Equations For Static Perfect Fluids With 
Spherical Symmetry,” J. Math. Phys. 28 (1987) 2949. 

[9] S. Rahman and M. Visser, “Spacetime geometry of static fluid spheres,” Class. Quant. Grav. 19 (2002) 935 
arXiv:gr-qc/0103065 . 

[10] K. Lake, “All static spherically symmetric perfect fluid solutions of Einstein’s Equations,” Phys. Rev. D 67 (2003) 104015 
arXiv:gr-qc/0209104 . 

[11] D. Martin and M. Visser, “Algorithmic construction of static perfect fluid spheres,” Phys. Rev. D 69 (2004) 104028 
arXiv:gr-qc/0306109 . 





